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Electromagnetic Potential Vectors 
and the Lagrangian of a Charged Particle 


John V. Shebalin 1 
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Abstract 

Maxwell’s equations can be shown to imply the existence of two independent three-dimensional 
potential vectors. A comparison between the potential vectors and the electric and magnetic field 
vectors, using a spatial Fourier transformation, reveals six independent potential components but 
only four independent electromagnetic field components for each mode. Although the 
electromagnetic fields determined by Maxwell's equations give a complete description of all 
possible classical electromagnetic phenomena, potential vectors contain more information and 
allow for a description of such quantum mechanical phenomena as the Aharonov-Bohm effect. A 
new result is that a charged particle Lagrangian written in terms of potential vectors automatically 
contains a ‘spontaneous symmetry breaking’ potential. 
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Introduction 


It is well known that the electric field E and the magnetic field H can be determined through 
the four-potential (<p,A). However, the scalar potential <p is merely an ‘auxiliary’ variable since, as 
will be shown here, it is only the three-vector A which is dynamically important. In transforming 
to a potential formulation, the first-order Maxwell's equations for E and H are replaced by second- 
order equations for the vector potential A. These second-order equations for A can be expressed 
as first-order equations for A and C, where C = 9,A (C and A must be treated as independent 
vectors, unless A is known for all time a priori ). Thus, a formulation of classical electrodynamics 
in terms E and H (the electromagnetic fields) is embedded in a formulation in terms of A and C 
(the ‘potential vectors’), which serves to describe electrodynamics in both classical and quantum 
mechanics. 

Maxwell's equations contain dynamic ( i.e ., time-evolution) equations only for the curls of 
E and H, while their divergences are prescribed. The equations for the potential vectors A and C, 
on the other hand, dynamically determine both their curls and divergences. If these systems of 
vector fields and partial differential equations are spatially Fourier transformed, then these 
preceding statements can be given in terms of the Fourier coefficients associated with the various 
spatial modes (which are identified by their wave vector k). Maxwell's equations dynamically 
determine only the transverse parts of the electromagnetic fields {i.e., kxE(k,t) and kxH(k,t)), 
while the longitudinal parts {i.e., k-E(k,t) and k H(k,t)) are prescribed. The potential vector 

equations, on the other hand, dynamically determine all the components associated with a given 
mode: kxA(k,t), kxC(k,t), k-A(k,t), and k-C(k,t). Maxwell's equations thus dynamically 

determine only four components for each mode (the transverse ones) which is sufficient 
classically, while the potential vector equations determine six components (both longitudinal and 
transverse) which is needed quantum mechanically. 

A direct, algebraic transformation from the set of modal vectors A(k,t) and C(k,t) to the set 
E(k,t) and H(k,t) will be given, and will show that the transformation is, in fact, a projection from 
(for each mode) a six dimensional space onto a four dimensional space. Again, a knowledge of 
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A(k,t) and C(k,t) for all modes gives a description of electromagnetic phenomena which 
encompasses both classical and quantum mechanics, while a knowledge of E(k,t) and H(k,t) for 
all modes gives a description which is restricted to classical mechanics alone. For example, an 
understanding of the Aharonov-Bohm effect requires that we recognize that A has a direct 
influence on physical processes; in classical electrodynamics, A is a convenient mathematical 
artifice, a potential which enters into physics only through its derivatives. 

In this paper, the preceding concepts will be presented in detail, utilizing spatial Fourier 
expansions. (Temporal Fourier expansions are not utilized because the coupling between the 
electromagnetic field and matter is, in general, non-linear, which precludes assigning a specific 
frequency to each spatial mode.) In the end, one important result will become clear: the 
Lagrangian of any system containing electromagnetically interacting matter must contain a term 
linera and a term quadratic in the particle. Since these charge density coefficients are themselves 
quadratic functions of quantum mechanical wave functions, the Lagrangian of any electrically 
charged particle must contain a term which is quadratic and a term which is quartic in the wave 
function associated with that particle. These naturally occurring terms are similar to the ad hoc 
‘spontaneous symmetry breaking’ potentials which have been introduced into current theories of 
elementary particles. 


Potential Vectors 

To begin the discussion here, consider Maxwell's equations: 

a) V D = p b) VxE = -3 t B 

c) V B =0 d) VxH =3 t D +j 

D = e 0 E, B = |i 0 H 

Here, p is the total charge density and j is the total electric current density. Although (1) is written 
in SI units, it will be more expedient to shift to ‘natural’ units where the electric permittivity is 
unity: £ c =l, as is the magnetic permeability: (i 0 =l. Maxwell’s equations then become: 
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a) V-E = p 
c) V H = 0 


b) VxE = -3 t H 
d) VxH = 9 t E +j 


( 2 ) 


Equation (2b) states that V H= constant and (2c) states that this constant is zero. Equation 
(2a) can be viewed as equating charge density p to V-E. The rest of E is dynamically determined 
by (2d); taking the divergence of (2d) and using (2a) we obtain the equation of continuity (or time- 
evolution equation for V-E): 


dtP+V-j=0 (3) 

The equations (2a) and (2d) ensure the existence of the equation of continuity, which has an 
obvious physical meaning; alternatively, (3) can be viewed as merely mathematically defining 
future values of V-E. 

In a similar manner, taking the curls of (2b) and (2d) can produce time-evolution equations 
for VxH and VxE, respectively. Presuming the values of E and H at infinity are known, then the 
Helmholz theorem [1] states that the six component equations found in (2b) and (2d) serve as 
independent, first-order, time-evolution equations which exactly determine the six components of 
E and H, as is mathematically required [2], 

The fields E and H can be represented in terms of a four-potential (<p. A): 

a) E = -Vcp -8 t A b) H = VxA c)a t <p+V-A=0 (4) 

These equations are well known: (4a) and (4b) satisfy (2b) and (2c) identically, while (4c) is the 
Lorentz condition, which arises because the divergence of A needs to be defined (the curl of A is 
defined by (4b)). However, the need for relativistic invariance requires the left-hand-side of (4c) 
be defined, rather than just V-A. The right-hand-side of (4c) can initially be an arbitrary (well- 
behaved) scalar function c, rather than zero; however, a gauge transformation of (cp. A)— K(p + 
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9,\|/, A-V\|/) can return the divergence condition to the form it has in (4c), as long as \\f satisfies the 
inhomogeneous wave equation 3 t 2 \j/ -V 2 \|/=a. Thus, the Lorentz condition can always be regained 
from an arbitrary choice of gauge. 

Using (3), and (4), Maxwell's equations (2a) and (2d) become: 


a) V 2 cp -d?<p = -p b) V 2 A -3? A = -j 


( 5 ) 


Here it would appear that the six (five non-trivial) independent evolution equations of (2) have 
been replaced by at most four equations; at least this is what is implied by stating that E and H are 
determined by cp and A. 

The source of the apparent over-determination is, of course, that Maxwell's equations (2) 
are first-order equations, while the equations (5) are second-order. Before changing (5) to a set of 
first-order equations, let us use (4c) and (5a) to obtain the following: 

V 2 cp = -p-V3 t A ( 6 ) 

Now, defining C = 3 t A, equations (5b) and (6) can be written as: 

a) V 2 cp = -p -V • C b)3 t C = V 2 A+j c)d t A = C (7) 


Thus we arrive at the realization that E and H can be found by determining A and C, and 
not just A (and <p) alone. The six independent evolution equations in (7b) and (7c) contain (2b) 
and (2d), and (7a) simply determines tp at each instant that p and C are known (p is presumably 
determined by using the equations of motion of whatever charged matter is present). If A is a 
known vector function, then C is determined through (7c) and only in that special case can E and 
H be found using A alone; otherwise, both A and C are needed for completeness. To continue 
this discussion, let us bring in Fourier representations at this point. 
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Fourier Representation 


The electromagnetic field and vector potential formulations of electromagnetism can easily 
be compared in terms of spatial Fourier expansions (here, the physical domain is assumed to be a 
finite-sized box): 


E(x) = (27t)' 3/2 X E(k)e ik x , etc. 

|k|<k™, (8) 

In (8), the number of modes for !kl< k mai is N, and the argument of E denotes whether we are in 
physical space or in Fourier space; in either case, time is omitted from the argument for brevity. It 
is presumed that the limit km ax — > °° may be taken and also that summation over discrete k can be 
transformed into a continuous integral over dk by expanding the size of the physical domain out to 
infinity. (Fourier expansions in time are not utilized since the system in which the electromagnetic 
fields arise must, in general be assumed to be non-linear, which does not allow for simple 
dispersion relations, i.e., a single frequency to be assigned to each mode k.) 

The various differential x-space relations can easily be transformed to the Fourier domain 
by the substitution V— >ik. Thus, the fields E and H, as given in (4), and A, C, and cp, as given in 
(7), can be algebraically related through their Fourier coefficients: 

a) E(k) = -ik(p(k) -C(k) b) H(k) = ikxA(k) 

c) <p(k) = k' 2 [p(k) +ik-C(k)] ( 9 ) 

In particular,each Fourier coefficient of the magnetic field H(k) is clearly due only to the transverse 
part of A(k); also, it is clear that k-A(k) does not play any role in determining the electromagnetic 

field components E(k) and H(k). 

The scalar function k-C(k), on the other hand, helps determine <p(k), as (9c) shows. 
Placing (9c) into (9a) yields: 
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a) E(k) = -ikk' 2 p(k) -(i -k 2 kk) C(k) 


b) H(k) = ikxA(k) 


00 ) 


where I is the unit dyadic. In this set of equations, it is clear that the transverse parts of E(k) and 
H(k) are determined only by the transverse parts of C(k) and A(k), respectively, and that ik-E(k) 
= p(k) and ik-H(k) = 0, as required. Apparently, the six components of A(k) and C(k) project 
onto only four components of E(k) and H(k). 

The dynamic Maxwell s equations (2d) and (2b), and vector potential equations (7b) and 
(7c), are, respectively: 

a) dE(k)/dt = ikxH(k) -j(k) b) dH(k)/dt = -ikxE(k) (n) 

a) dC(k)/dt = -k 2 A(k) +j(k) b) dA(k)/dt = C(k) (12) 

Taking a dot product of (10a) and (1 la) with ik and combining the results yields the modal form of 
the continuity equation (3): 


dp(k)/dt + ikj(k) = 0 


(13) 


This is not so much an evolution equation as an identity which must be satisfied by the dynamic 
equations of both the electromagnetic fields and charged matter. 

The Lorentz condition (4c), in terms of spatial Fourier expansions, appears as: 


dtp(k)/dt + ik-A(k) =0 (14) 

Using (9c) and (14) yields (13), as it should. Also, a gauge transformation has the modal form: 
{ A(k), C(k)} — { A(k) - ik vj/(k ), C(k) - ikd\|/(k)/dt}, where \j/(k) satisfies the modal wave 
equation d 2 \j/(k)/dt 2 + k 2 y(k) = 0. 
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In (12), the occurrence of two additional dynamic components in A(k) and C(k) , 
compared with E(k) and H(k), indicates that a vector potential formulation may have more explicit 
physical information than a formulation in terms of electromagnetic fields alone. In particular, (10) 
shows that neither V A or V C determine E or H. However, the longitudinal parts of A and C 
and not just their transverse parts, may have observable effects. This implication has, in fact, been 
realized for the longitudinal part of A in the Aharonov-Bohm effect [3]. Before discussing this and 
other matters, however, let us consider the completeness and incompleteness of the various 
dynamical representations. 

Complete and Incomplete Dynamical Representations 

So far, it has been shown that underlying Maxwell's equations there is a potential vector 
formulation of electromagnetics. In particular, it was shown that in addition to the vector potential 
A, we must have another independent vector C to complete the formulation. This additional vector 
would, in fact, be equal to the partial time derivative of the original vector potential, if that vector 
were already known. Otherwise, it must be treated as completely independent, similar to the 
manner in which x and v are treated in the Lagrangian formulation of dynamics, or x and p in the 
Hamiltonian formulation of dynamics. 

The analogy to the canonically conjugate variables of particle mechanics can be made more 
precise in the Fourier representation. Rather than the position x ; and momentum p ; of the i* 
particle serving as canonically conjugate variables, the coefficients of the electromagnetic field E(k) 
and H(k) of the k th mode are the canonically conjugate variables. In the vector potential 
formulation the canonically conjugate variables are A(k) and C(k). (In either of these 
formulations, it is the source vector j(k) which ultimately give rise to the coupled vector fields, and 
also ties them into material motion.) 

The analogy with particle mechanics can be extended further: in the same way that the set 
of positions and momenta {Xj, P; I i=l,...,N) defines a 6N-dimensional phase space, the set of 
conjugate variables (A(k), C(k) I lkl<k m „) on N modes also defines a 6N-dimensional phase 
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space. Here k and -k are counted as separate modes even though, since A(x) and C(x) are real, 
A(k)=A*(-k) and C(k)=C*(-k). The complex vectors A(k) and C(k) are composed of 6 
components each ( i.e ., A(k)= A R (k)+iAi(k), where Arj are real) and the pair (A(k), C(k)) thus 
contain twelve ‘degrees-of-freedom’; however, these are the same as are found in the pair (A(-k), 
C(-k)). Then only half of the N modes, such that lkl< are independent and 12xN/2=6N. 

The phase space associated with the set {E(k), H(k); Ikl^k^J of N modes has, however, 
a dimension of only 4N because k-E(k) and k-H(k) are prescribed. Transformation (10), which 
takes (A(k), C(k); lkl<k m „}^{E(k), H(k); lkl<k m3I ) is actually a projection operation and is not 
invertible. Although the set of vectors (E(k), H(k); Ikl^k^J may be sufficient classically, it is 
not, as has been mentioned, sufficient quantum mechanically. The classical phenomena described 
by E and H thus fall into only a subset of all possible electromagnetic phenomena. 

The Aharonov-Bohm Effect 

Consider an infinitely long wire with a steady electric current ‘I’ coursing down it. In 
order to determine the magnetic field, we solve Maxwell's equation (2d) in a time-independent 
form, VxH = j, where j is zero outside the wire. The solution is elementary and well-known: in a 
cylindrical coordinate system (r,<j>,z), whose z-axis runs down the center of the current carrying 
wire, outside of the wire the only non-zero component of H is H < ,=I(27tr) 1 . 

Consider, instead of a long current-carrying wire, a long, straight, and narrow solenoid. 
This solenoid will completely contain a magnetic field of total flux d>; the vector potential A will 
now satisfy the equation VxA = H, where H is zero outside the solenoid. Since this situation is 
mathematically the same as the current-carrying wire problem of the previous paragraph, we can 
immediately give the solution outside the solenoid: the only non-zero component of A is 
A < ,=d>(2ttr) 1 . 

Classically, an electron beam passing by the long solenoid would feel no electromagnetic 
influence from the solenoid, since its motion is only affected by E and H, which are zero outside 
the solenoid. Quantum mechanically, however, the electron's wave function is affected: 
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Depending on which path we trace the electron beam's motion along, the phase 0 of the wave 
function associated with the beam has an additional value of (here Planck's constant has a value of 
unity) 


A0 = e 



Adi 


(15) 


compared to what it would have if A (due to the flux d>) were absent [3]. 

In fact, a solenoid placed between the beams of a double slit electron diffraction experiment 
will cause a different phase shift to occur for the wave function of each beam. The difference in 
phase shift for the wave functions of the beams, one through the right-hand slit and the other 
through the left hand slit, with the solenoid between them, will be given by the application of 
Stokes theorem to (15) around a closed path: 


A ©right ' A ©left — 6 



VxA-ds = ed> 


(16) 


Thus the phase difference, which will cause the diffraction pattern to shift as the solenoid's field is 
varied, depends only on the total magnetic flux between the beams. This effect has been observed 
in a number of careful experiments [4], 

The vector potential outside of the solenoid due to the flux d> satisfies VxA=0. Thus 
kxA(k) = 0 and it is only the part of A(k) along k which can cause the effect, i.e., it is k-A(k) 
which is important here. Since this is precisely what is projected out of the electromagnetic fields, 
an effect such as that of Aharonov and Bohm can be overlooked if only E and H are assumed to 
have physically observable effects. 
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Charged Particle Lagrangians 


At this point, let us consider what novel effects are uncovered when the Lagrangian for a 
charged particle is expressed in terms of potential vectors. The Lagrangian density for any panicle 
of classical or quantum mechanical electric charge e will contain a part A e related to the 
electromagnetic field: 


A e = -eptp + ej-A + L(e 2 -H 2 ) 


(17) 


(Note: The symbols (p, j) will now be used to signify the particle four-current density; the electric 
four-current density will now have components (ep, ej).) 

The Lagrangian itself is defined as a volume integral of a Lagrangian density: 


L(t) = 


I 


A(x,t) d 3 x 


(18) 


Since we have defined the system volume to be a periodic box with an associated Fourier 
expansion (8), the Lagrangian could equally well be written in terms of Fourier coefficients: 

L(t) = £ A(k,t) 

k (19) 

where the k-space Lagrangian density is (surpressing explicit time dependence for brevity): 

A(k) = A m (k) - ep*(k)q>(k) + ej*(k)A(k) + l(| E(k)| 2 -| H(k)| 2 ) 

In (20), A m contains those terms which pertain only to the matter in the system; the specific form 
of A m depends on whether a classical or quantum mechanical system is under condsideration. The 
rest of the Lagrangian pertains equally to both classical or quantum mechanical systems. Notice 
that if we set (E, H) equal to zero in (20), then the Lagrangian density is A(k) = A m (k) - ej^A^, 



i.e., it still depends on the four-potenial. Thus, even though (E, H) are zero, the particle is not 
free (for which the four-potential must also be zero), and effects such as that of Aharonov and 
Bohm are possible. 

If the expressions for <p(k), E(k), and H(k) ((9c), (10a), and (10b), respectively) are 
placed into (20), the result may be written as: 

A(k) = A m (k) + A r (k) - A p (k) (21) 


where ‘radiation’ part \ and the ‘potential’ part Ap of A c are: 


a) A r (k) = ^-(|C±(k)| 2 - k 2 |A x (k)| 2 ) 


b) A p (k) = p*(k)(i(k) - ej*(k)-A(k) +|-k- 2 |p(k )| 2 


( 22 ) 


Here C x and A x are the transverse parts of C and A, respectively. The coefficient |i(k) is: 

p(k) = iek' 2 k-C(k) (23) 


Here we see that k-C(k), though it does not appear in the electromagnetic fields, does appear in 
the Lagrangian density. 

In the Lagrangian density (22), the ‘potential’ A p (k) contains a term which is linear and one 
which is quadratic in p(k): p(k)p*(k)+Q(k), where Q(k) = ie 2 k- 2 |p(k)| 2 . What does Q(k) look 
like in x-space? If we define a(k) = p(k)/k for k>0 and a(k) = 0 for k=0, then we have: 


y ^|o(k)| 2 = I d 3 x^[a(x)] 2 

k J 


( 24 ) 


Here a(x) is the inverse Fourier transform of o(k), etc. Thus Q(x) = ^e 2 [a(x)] 2 . 


Although the potential (22) is written in k-space, it has a striking resemblance to ad hoc 
‘spontaneous symmetry breaking’ potentials found in current theories of particle physics [6]. (In 
these theories, p(x) is the modulus of the quantum mechanical wave function squared.) The 
naturally occuring term we have uncovered here would appear to deserve further study. 

Conclusion 

In this paper, four major points have been discussed. First, there is the potential vector 
formulation of electromagnetic fields. It was seen that the divergence of the vector potential is 
uniquely defined, up to a gauge transformation. The potential vector formulation contains 
Maxwell s equations explicitly, while it might be said that Maxwell's equations, Lorentz 
invariance, and gauge invariance, all taken together, implicitly contain the potential vector 
formulation. 

The second point was that of dynamical completeness. For a given number N of Fourier 
modes, it was clear that the electromagnetic fields associated with these N modes span a phase 
space which is only 2/3rds the dimension of that spanned by the potential vectors for the same 
number of modes. Thus, the electromagnetic fields offer only an incomplete dynamical 
representation, as compared to the potential vectors. 

The third point was to review how this incompleteness could miss an observable effect, 
albeit a quantum mechanical one. This observable effect is that of Aharonov and Bohm, where the 
longitudinal part of A proved important. This clearly demonstrates how fruitful it is to have and to 
use a conceptually complete dynamical representation of electromagnetic phenomena. 

The fourth and final point was to recognize the existence of a naturally occuring 
‘spontaneous symmetry breaking’ term in a charged particle Lagrangian. This term contains the 
longitudinal part of C, indicating that it is also dynamically important. The detailed effect of this 
term on quantum mechanical descriptions of charged particle dynamics is a separate issue from the 
topics covered herein, though clearly an intriguing one. 
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